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Integrating out the Dirac sea: Effective field theory approach 
to exactly solvable four-fermion models 
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We use 1+1 dimensional large N Gross-Neveu models as a laboratory to derive microscopically ef- 
fective bagrangians for positive energy fermions only. When applied to baryons, the Euler-bagrange 
equation for these effective theories assumes the form of a non-linear Dirac equation. Its solution 
reproduces the full semi-classical results including the Dirac sea to any desired accuracy. Dynamical 
effects from the Dirac sea are encoded in higher order derivative terms and multi-fermion interac- 
tions with perturbatively calculable, finite coefficients. Characteristic differences between models 
with discrete and continuous chiral symmetry are observed and clarified. 
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I. INTRODUCTION 

Exactly solvable models in quantum field theory which 
bear any resemblance to the real world are extremely 
rare. One example is provided by the Gross-Neveu 
model family pj, 1+1 dimensional four-fermion interac- 
tion models with N flavors and Lagrangian 

C = $ (i0 - m ) V + y [W) 2 + A(# 75 </>) 2 ] , (1) 



where A = or 



and flavor indices are suppressed 
(iptp — V'feV'fe etc.). Depending on the value of 

A, these models feature discrete (if) — > 75^, A = 0) 
or continuous [tp — > exp( 10:75 )V->, A = 1) chiral sym- 
metry, possibly broken by the bare mass term ~ tuq. 
To avoid confusion we shall refer to the first model as 
Gross-Neveu (GN2), the 2nd one as Nambu-Jona-Lasinio 
(NJL2) model. From the point of view of strong inter- 
action physics, they are most useful in the 't Hooft limit 
(N — > 00, Ng 2 = const.) to which we stick in the follow- 
ing. Rather than repeating all the well-known attractive 
features of these models, we refer the reader to some per- 
tinent review articles 0, HI, H| and the references therein. 

The Lagrangian |T]) is known to possess multi-fermion 
bound states analogous to baryons or baryonium states in 
hadron physics. They were first constructed by Dashen, 
Hasslacher and Neveu (DHN) in the massless GNo model 
H and by Shei in the massless NJL2 model [6f. The 
semi-classical method developed by these authors may 
be rephrased equivalently as a Dirac-Hartree-Fock calcu- 
lation @, [!]• Here, one solves the first-quantized, time- 
independent Dirac equation 



(-75^ + 7 °S + rfP) = E a i> a 



(2) 



for single particle orbits with label a subject to self- 
consistency conditions for scalar and pseudoscalar mean 



fields, 



S = mo - Ng 2 ^2 ipaiia , 

a 

occ 

P = -\Ng 2 ^2iP a i l5 i> a . 



(3) 
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A major challenge arises from the fact that the Dirac sea 
must be included in the sum over occupied states in ((3|) . 
This problem was solved in Refs. [5|,[6| by inverse scatter- 
ing methods together with a careful subtraction of bound 
state and vacuum energies. At about the same time the 
proposal was made to solve these models "classically", 
neglecting the Dirac sea and including only the discrete 
valence level [1, [l(| ■ This reduces the full problem to 
a non-linear Dirac equation, which has been solved in 
closed analytical form for both the GN2 and NJL2 mod- 
els. It appeared that the classical fermionic solution was 
useful for the GN 2 model but completely failed in the 
NJL2 case @ for poorly understood reasons. The Dirac 
sea also caused considerable difficulties in attempts to 
base nuclear physics on field theoretic models like the 
(<j,Cl>) model [Hfl . so that mean field calculations were 
typically done without the sea. 

In the present work we reconsider the role of the Dirac 
sea in multi-fermion bound states of Gross-Neveu models. 
From a general field theoretic point of view, it should be 
possible to "integrate out" the Dirac sea and derive an ef- 
fective Lagrangian for (positive energy) valence fermions 
only. In the large N limit, it would then indeed be suffi- 
cient to solve the classical Euler-Lagrange equation of 
this effective theory for fermions. From such a point 
of view the above mentioned "no-sea" calculations may 
be re-interpreted as follows: The authors assumed that 
the effective Lagrangian is the same as the original La- 
grangian (P), except that the bare parameters mo and g 2 
are replaced by an effective mass and coupling constant. 
In the course of this paper we will confront this implicit 
assumption with our results for an effective action de- 
rived from the underlying field theory. 

Let us mention that although the solvable models used 
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here are restricted to 1+1 dimensions, the technique ap- 
plied in deriving the effective Lagrangian is not. It might 
be useful also in higher dimensional (mean field type) the- 
ories. The advantage of developing the methods in the 
context of Gross-Neveu models is the fact that the ex- 
act bound states are known analytically. This enables 
us to test our effective action quantitatively and make 
sure that the expansion is consistent to a given order in 
some small parameter. Since the derivation of £ e ff is not 
completely straightforward due to the necessity of resum- 
mations, this has turned out to be quite helpful indeed. 

This paper is organized as follows. In Sect. II we derive 
the effective Lagrangian for the GN2 model in the three 
lowest orders of a systematic expansion valid for small 
filling fraction and/or large bare fermion mass and check 
it against the DHN baryon. In Sect. Ill we take up the 
NJL2 model. Due to additional complications, we will 
content ourselves with the leading order effective action 
here and test it against the Shei bound state. We end 
with a short summary and conclusions in Sect. IV. 
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FIG. 1: Dyson equation for fermion propagator in HF ap- 
proximation. Dashed line: free propagator, solid line: dressed 
propagator. 
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FIG. 2: Dynamically generated fermion mass in the vacuum. 
A typical "cactus" diagram produced by an iterative solution 
of the HF equation in Fig. 1 is illustrated. 



II. GROSS-NEVEU MODEL (GN 2 ) 

Our aim is to account for the effects of the Dirac sea by 
means of an effective theory of positive energy fermions 
only. Since we do not use the path integral but work 
canonically, we first have to explain what we mean by 
"integrating out the Dirac sea" . We start with the mass- 
less GN2 model described by the bare Lagrangian 



(4) 



and focus on the DHN baryons. They are character- 
ized by a discrete positive energy level occupied with 
n < N fermions, all negative energy levels being com- 
pletely filled. The presence of extra fermions polarizes 
the Dirac sea, distorting the single particle wave func- 
tions. Ideally, one would like to fully integrate out the 
Dirac sea. We have not been able to do this. Anyway, 
the full effective theory will invariably be non-local and 
therefore of less practical use. We therefore look for an 
expansion parameter which would enable us to derive an 
"almost local" effective action, consisting of a polynomial 
in the chiral condensate tpip and its first few derivatives. 
From the DHN baryon we know that for small filling 
fraction v = n/N of the valence level, the self-consistent 
scalar potential S(x) differs from the physical (vacuum) 
fermion mass to by a weak and slowly varying potential 
only. This suggests to use v as expansion parameter and 
to restrict oneself to small filling fraction where the HF 
potential becomes soft. 

The following derivation of the effective action is tai- 
lored to the HF approach to which we now turn. In 
the vacuum, a fermion mass is generated dynamically as 
shown graphically in Fig. 1. The self-consistent mass 
can be thought of as the sum over all one-particle- 
irreducible (1PI) cactus type diagrams, see Fig. 2. The 



(self-consistent) tadpole diagram receives contributions 
from all negative energy occupied states, leading to the 
gap equation 



to 1 



Ng 2 A 







(5) 



(with A/2 as UV cutoff Q). Clearly, since the physical 
fermion mass is a pure manifestation of the Dirac sea, it 
has to be put in by hand into an effective theory as an 
effective mass term 



,(i) 

-off 



(6) 



(From now on the superscript on C e s refers to the number 
of loops of the self-energy diagram from which it was 
derived. Due to resummations, it may contain higher 
loop effects as well.) 

Let us now turn to the problem of finite fermion num- 
ber. The HF approach still has the same basic structure 
as in Fig. 1, except that the fermion self-energy is in 
general x-dependent. The propagator gets an extra con- 
tribution from the positive energy valence states. We 
denote the vacuum and valence particle contributions by 
"— " and "+" , respectively. For the free massive propaga- 
tor for instance, the corresponding decomposition can be 
inferred from the well known result for finite temperature 
and chemical potential [HI, [l3j] , 



sip) 

S-(P) 

s+(p) 



S-(p) + S+(j>), 
i 



j> — to + ie ' 
= — 2n5(p 2 — m 2 )(fi + to) 
x [6{- Po )9(E f - E p ) + 6( Po )6(-E f - E p 



(7) 



Here, S- is just the free Feynman propagator, Ef the 
Fermi energy or chemical potential. The one-loop self- 
energy (i.e., the tadpole) then naturally splits up into the 



FIG. 3: Decomposition of the fermion self-energy into a) neg- 
ative and b) positive energy contributions. 




FIG. 4: Topologically distinct diagrams contributing to the 
fermion self-energy in HF approximation from two up to five 
loops. They are needed to construct the effective action for 
the GN2 model in NNLO, see the main text. 

two pieces shown in Fig. 3. In the effective theory with 
the Dirac sea integrated out, diagram a) is accounted 
for by the mass term ([5]). Diagram b) corresponds to 
the lowest order "no-sea" HF calculation with the origi- 
nal four- fermion interaction from Eq. (fj|. This splitting 
of the tadpole into +/— pieces is also the key for the 
following systematic procedure. For a given number of 
loops, draw all topologically distinct 1PI cactus diagrams 
as shown in Fig. 4 up to five loops. Split each loop into 
a sum of " + " and " — " contributions. This generates 
2 n distinct, "labeled" diagrams out of each n-loop dia- 
gram in Fig. 4. Then analyze each labeled diagram with 
respect to the question whether it is accounted for by 
the effective action with the Dirac sea integrated out. If 
it is not, add a new term to £ c g which generates the 
corresponding fermion self-energy. 

Although this scheme is straightforward and construc- 
tive, it is still too naive. We will see that it would 
be meaningless to terminate the procedure at any finite 
number of loops. Remember that the mass term (j6j) al- 
ready includes an infinite number of diagrams due to the 
self-consistency condition (see Fig. 2). Otherwise, di- 
mensional transmutation, a truly non-perturbative phe- 
nomenon, could not occur. We can only make sense out 
of the n-loop contribution to the effective Lagrangian if 
we perform appropriate resummations. This turns out to 
be a crucial element of the whole procedure and the way 
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FIG. 5: Labeled diagrams selected from the diagrams in 
Fig. 4, which are the seed for new terms in the effective La- 
grangian. See main text for a detailed discussion. 
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FIG. 6: Resummation of scalar vacuum polarization graphs 
into a momentum dependent effective coupling. The infinite 
sum of " — " bubbles may be thought of as generating the a 
propagator. 

in which the bare coupling constant g 2 gets traded for a 
finite effective coupling constant. 

To explain how this works, let us go through the first 
few terms of the loop expansion in some detail. Guided 
by the DHN baryon at small filling, we shall treat 'tpip 
and k 2 (or □) as being of the same order e and derive 
systematically all terms up to 0(e 4 ). 

At two-loop order, the only diagram (out of 4 labeled 
diagrams) which induces a new term in £ e g- is the one 
shown in Fig. 5a. This corresponds to a vertex correc- 
tion to the " + " tadpole in Fig. 3b. Technically, it in- 
volves the scalar vacuum polarization loop (labeled " — " 
in Fig. 5a), where scalar means that the vertices are just 
1. We evaluate it by standard Feynman rules and ex- 
pand the result in powers of the 4-momentum k flowing 
through the graph. The result for the sum of the " + " 
tadpole and Fig. 5a then gives the following self energy 
contribution, 

SS(k) = - 5 2 (^)fe{l + ^-fln--l 

( 7T \ TO 

Here, the scalar condensate in momentum space, (Tptpjk, 
is defined so as to include only " + " terms. Whereas the 
fc 2 -dependent correction terms are suppressed (remember 
that k 2 ~ e), the k = contribution ~ ln(A/m) needs 
a special treatment. According to the gap equation ([5]), 
this term is actually of 0(1) so that the series expansion 
starting with Eq. ([5]) is meaningless as it stands. We 
therefore sum the whole series of graphs shown in Fig. 6 
into a geometrical series, or, equivalently, replace the bare 
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coupling constant by the effective, momentum dependent 
coupling 



SoffW 



9 



Ng 2 



In — — 1 



k 2 
12 m 2 



+ 



(k 2 ) 2 
120m 4 



(9) 



Due to the gap equation, the 1 in the denominator is now 
cancelled against (Ng 2 /ir) ln(A/m), and the bare cou- 
pling constant g 2 drops out of this expression. Expand- 
ing the ^-dependent terms again [to 0(fc 4 )], we arrive at 
the following momentum dependent effective coupling, 



N 



ll(k 



2\2 



12m 2 



720m 4 



(10) 



It is now easy to write down a two-loop effective La- 
grangian which would give just these correction terms 
(correct at two-loop level, but containing higher order 
terms due to resummation), 



/•(2) _ 



2N yvV > 
IItt 



1440iVm 4 



24iVm 2 
(□ 2 ^)(^) 



(□^)(^) 



(11) 



Notice that the bare coupling constant g 2 in the original 
(ipip) 2 term has been replaced by ir/N, and new four- 
fermion interaction terms containing derivatives of rptp 
have been generated. Higher order terms could easily 
be obtained from the full fc-dependence of the vacuum 
polarization graph if desired. The result is manifestly 
non-perturbative, as it does not contain g 2 , and finite 
owing to the use of the gap equation. 

It may be worthwhile to pause here and comment on 
the value of the leading order effective coupling, g 2 s = 
tt/N [sec Eq. (fTU|) ] . This quantity already appeared in the 
original paper by Gross and Neveu [l[ as fermion-fermion 
scattering amplitude at zero momentum. These authors 
also discuss the scalar a meson, pointing out that the 
square of the fermion-antifermion-u coupling constant is 
given by g 2 F p = Airm 2 /N. One can then understand the 
effective coupling in a similar way as in the old Fermi 
theory of weak interactions, namely as a product of two 
coupling constants and a heavy boson propagator (M a = 
2m), 



9aFF 



1 

M 2 ' 



jFF 



71 

N' 



(12) 



Keeping the momentum dependent terms approximately 
accounts for the finite range of the cr-exchange. 

We now turn to the three-loop graphs. The two topo- 
logically distinct graphs in Fig. 4b give rise to 2 x 2 3 = 16 
labeled subgraphs. By inspection we find that only the 
graph shown in Fig. 5b is the seed for a new term in the 
effective action, all other diagrams being generated by 
lower order terms or self-consistency in the " + " sector. 
Evaluating the fermion loop with three scalar vertices 



and finite momenta k\, k 2 from the two " + " loops, we 
find the contribution to the self-energy 



(5S(fc) 



1 + 



dfci dk 2 

k\ + k\ 



(27r)<5(fc 
: + kik 2 



(Ng 



2\3 



6m 2 



ki - .._ , 

27rm 

(#)ti(#)t ; 
N 2 



(13) 



Here, it is sufficient to go to 0(k 2 ). Resumming bub- 
bles by replacing each of the three couplings Ng 2 at the 
vertices by effective ones analogously to Fig. 6, 

(Ng 2 f -> Ng^{k 1 )Ng^(k 2 )Ng^(k 1 + k 2 ), (14) 

we find to 0(fc 2 ) 

6Z(k) = - I ^^l(2n)6(k-h-k 2 )f- (15) 
J 2tt 2tt 2m 

fc 2 +fc|+fclfc 2 \ (^V)fe 1 (^'i/ , )fc2 



3m 2 



N 2 



The effective Lagrangian that yields this self-energy in- 
cludes the following six-fermion interactions, 



C 



(3) _ 



6mN : 



12m 3 iV 2 



(□Vx/OW^) 2 - (16) 



At four-loop order, there is again only a single graph 
(out of 5 x 2 4 = 80 labeled diagrams, see Fig. 4c) which 
generates a new term in the effective action. It is shown 
in Fig. 5c. The central loop labeled " — " with four scalar 
vertices is easy to compute, since we only need its value at 
k = 0. Resumming the bubbles by substituting Ng 2 — > ir 
[to O(fc )], an eight- fermion interaction, 



24m 2 N 3 



(17) 



is induced. This is not yet the whole story to 0(e 4 ) 
though. Out of the 11 x 2 5 = 352 labeled five-loop di- 
agrams derived from Fig. 4d, the two shown in Fig. 5d 
give rise to a further eight-fermion interaction term of 
the same order as (117[) . Applying the by now familiar 
resummation it is given by 



,(5) = 7T 

"" eff 8m 2 Af 3 



(18) 



The different origin of the contributions (jTTJ) and JT5|> 
becomes clearer if one interprets the scalar bubble sum 
as cr-meson propagator. While the four-loop contribution 
(jTTJ) corresponds to a aa — > aa point-like interaction due 
to a heavy fermion loop, the five-loop term (fTS]) has a 
different topology. It arises from the process era — > a — > 
aa with an iterated effective three-cx vertex, see Fig. 7. 

This completes the discussion of all terms up to 0(e 4 ). 
We summarize by writing down the effective low-energy 
Lagrangian for the massless GN2 model valid to 0(e 4 ) 
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FIG. 7: Illustration of the difference between the eight- 
fermion interactions and in Eqs. (|I7|) and (|18fl . The 
wiggly lines are a mesons and represent the bubble sum shown 
in Fig. 6. Each external line ends with a scalar density ipip. 



By construction, the condensate (ipip) that appears here 
refers only to the positive energy sector. Using the matter 

part S+ of the Dirac propagator (p/ = 

denotes the Fermi momentum), we find 



M 

(i/ji/;) = N — In 

7T 



Pf 



Pf 



M 



M 2 



(23) 



where the Dirac sea has been integrated out, 



+ 



6mA 2 



llTT 



1440m 4 A 



(□ 2 ^)(#) 



12m 3 A 2 



6m 2 A 3 



OHO 4 - (19) 



By counting powers of tpip and □, we identify the first 
two terms as leading order (LO), the next two as next-to- 
leading order (NLO) and the remaining three as NNLO 
approximation. To LO, a mass term is generated and 
the coupling constant g 2 of the four-fermion interaction 
is replaced by the effective coupling constant tt/N. In 
higher orders, the Dirac sea manifests itself through mo- 
mentum dependent couplings and the appearance of six- 
and eight-fermion interactions, to the order we are work- 
ing. Here we are in the fortunate position of being able to 
evaluate all coefficients of the effective Lagrangian from 
the underlying field theory. Technically, the calculation 
of the coefficients only involves standard one-loop Feyn- 
man diagrams in the vacuum, without any reference to 
finite fermion density or baryons. Nevertheless, the La- 
grangian (p~9j) should be adequate to predict properties of 
baryons or baryonic matter by means of a purely classical 
calculation. 

We shall test this conjecture against known results for 
the full GN2 model in two different ways. First, we do 
a "no-sea" HF calculation for matter at low density, as- 
suming unbroken translational invariance. Second, we 
evaluate the baryon with small occupation of the valence 
level by solving a non-linear Dirac equation. 

If the condensate tyip) is assumed to be translationally 
invariant, only the non-derivative terms of the effective 
action enter, i.e., 



r' 



ip(i(j) — m)if) - 



2N 



W) 2 + 



6mN 2 



, ..,-,>■ • - ■ (20) 
bm/iV' 3 

The Euler-Lagrange equation allows us to identify the 
effective mass M via 



with 



M 



N 



2mN 2 



2tt 3 



3m 2 A 3 



(21) 



W) 3 - (22) 



A low density expansion in pf yields 



M = m — Pf 



2m 2m 2 



P} 



± ... 



(24) 



for the solution of Eqs. Q22I23|) . in perfect agreement with 
the Taylor expansion of the exact result [2|, 



M 



y m 2 — 2mp f . 



(25) 



As is well known, this particular HF solution is not a 
minimum of the action. If translational invariance is as- 
sumed to be unbroken, there is a first order phase transi- 
tion which requires a Maxwell construction, see Ref. 0. 
Nevertheless this unphysical solution can be used as a 
purely algebraic test of the effective action. We thus 
confirm that the effects of the Dirac sea are encoded in 
the effective action, both through the mass and coupling- 
constant of the four-fermion interaction and through in- 
duced many-fermion interaction terms. Obviously, the 
lower the density one is interested in, the smaller the 
number of terms that need to be kept in £ e g . 

A physically more relevant test case is the DHN 
baryon. We start with the Euler-Lagrange equation de- 
rived from the full effective Lagrangian (|19p and every- 
where replace ipifi by its ground state expectation value, 
this time determined by a single valence state. Let us 
define the ^-dependent condensate 



N 



(tpip) = -irvip ip 



(26) 



with v = n/N. Since s(x) is time independent, we get 

(27) 



N 

□ (^) = — s". 

7T 



Consequently, the Euler-Lagrange equation can then be 
cast into the form of the following non-linear Dirac equa- 
tion for the valence level, 



where the scalar potential S is given by 



(28) 



S = m + s + 



1 



1 



11 



jv 



12m 2 



2m 720m 4 
' [2*"s + (s') 2 ]+is 3 , (29) 



6m 3 



3m 2 
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and s in turn depends on ipo through Eq. 
Eqs. (|28I29|) have to be solved self-consistently subject 
to the normalization condition 



dxVwo = 1- 



(30) 



The mass of the baryon in the effective theory is then 
calculated as follows: Deduce the Hamiltonian density 
from C c s (|19j) in the standard way, 



"Heff 



dC 



dip 



4> - Asfl ■ 



(31) 



Since C e g is linear in ip, this amounts to dropping the 
term containing the time derivative of ip and reverting 
the overall sign of £ c ff- The baryon mass is now simply 
given by 



M B = / dxH eS - 



(32) 



To solve the "no-sea" Dirac-HF equation is not easy and 
would have to be done numerically in general. Since we 
know the exact answer for tjjo from the full GN2 model 
and are primarily interested in checking our effective ac- 
tion, we proceed differently. We use the following repre- 
sentation of the 7 matrices, 



7° = 



-01, 



7 =10-3, 



75 = 7°7 1 = 



-cr 2 . 



(33) 



The (positive energy) valence level of the DHN baryon 
then has the spinor wave function and energy [7J 



cosh 
1 



E Q = m\/l - y 2 , 



(34) 



shows that the terms retained as well as the resumma- 
tions done are consistent. 

In the massless GN2 model, the only regime where an 
almost local effective action for the baryon can be jus- 
tified is v <C 1, i.e., small valence filling fraction. Once 
we include a bare mass term, there is yet another handle 
to suppress Dirac sea effects and make the scalar poten- 
tial softer, namely by increasing the bare mass. This 
incites us to generalize the "no-sea" effective action to 
the massive GN2 model. The only difference to the pre- 
vious calculation comes from the modified gap equation 
which now reads UM 



Ng 



7T A 

= 7 + In — . 



While the diagrams used to derive the effective action 
are the same as above, the algebra slightly changes. We 
immediately jump to the final effective Lagrangian for 
the massive GN 2 model in NNLO, 



1 



7T 1 

24m 2 N (1 + 7) 2 



2iV(l+ 7 ) 



6mN 2 (1 + 7 ) 3 



W) 3 



llTT (1+ 67/H) /ri2 



(□VVOW) 



1440m 4 iV (1 + 7 ) 3 



+ 



(1 + 7 /4) 



6m 2 N 3 (I + 7) 5 



W) 4 



(39) 



with the definitions 



ymx ± — artanhy, 



V 



T- (35) 



We plug this ansatz into the non-linear Dirac equation 
(|28p. expand in the filling fraction v and check with 
MAPLE that the equation is indeed satisfied up to (in- 
cluding) 0(v & ). For the baryon mass (l32|) we obtain at 
this order 



M R = 



1 - 



24 



1920 



322560 J 



(36) 



We have truncated the series because higher order terms 
are not reliable, as they would require an improved effec- 
tive action. The exact DHN baryon mass is given by the 
compact formula 



Mb = nm- 



sin6* 



(37) 



with 8 as defined in Eq. (|35|) . If we expand this function 
in powers of v, the series indeed starts with (|36|) . This 
is obviously a very good test of our effective action and 



It reduces to Eq. (|T9|) in the chiral limit or, equivalently, 
at 7 = 0. In order to test our effective Lagrangian, we 
turn to the known exact baryons of the massive GN 2 
model [H [M [l|]. Defining s(x) as in Eq. (USD, the 
non-linear Dirac equation (|28p now contains the scalar 
potential 



S 



1 



1 



1 



(I+7) 1 12m 2 (1 + 7 ) 2 

1 1 2 11 (1 + 67/H) xv 



2m(l+7) 3 720m 4 (1 + 7 ) 3 



(40) 



1 (1+7/2) ro „„„ , ,_„ ai , 2 (1 + 7/4) 3 



6m 3 (I+7) 4 



[2 S ".s+( S ') 2 ] + 



3m 2 (I+7) 5 



The exact spinor wave function and the energy remain 
the same as in Eqs. (|3"4")l , but the relationship between 
the parameter y and the valence occupation fraction v 
changes to 



7 



tan( 



y = sin f 



(41) 



Eq. (|4T[) can be solved for 8 by means of a power series 
expansion in v for given 7. The exact baryon mass is 



7 



given by 



M R = 



2mN 



sin 8 + 7 artanh (sin 6)] 



(42) 



We have verified with MAPLE that the full wave func- 
tion ?/>o solves the non-linear Dirac equation derived from 
Eq. (j3"9"|) to an accuracy of 0(i/ 6 ), as in the massless case. 
The baryon mass calculated from the effective action is 
found analytically to be 

tt 2 v 2 7r 4 (l + 9 7 K 



Mb = nm - 



24(1 
■ 54 7 - 



- 7 ) 2 1920(1 
2257 2 )^ 6N 



. 7 )E 



322560(1 + 7 ) s 



(43) 



the generalization of Eq. (f3"6")> to arbitrary 7. To the given 
order this agrees once again with the result based upon 
the series expansion of the exact equations (|41l42j) . 

These results show that with our "no-sea" effective ac- 
tion, we have indeed achieved what we were aiming at. 
Depending on the number of terms one is willing to in- 
clude, one can systematically get an increasingly accu- 
rate value for the baryon mass and valence wave func- 
tion. Note that the expansion parameter changes from v 
in the chiral limit to 1//7 in the heavy quark limit. As a 
matter of fact, the truncated series (|43[) is a much better 
guide to the exact baryon mass in the full range of the 
parameters (y, 7) than could have been expected on the 
basis of our derivation. We find that the NNLO effective 
action yields the binding energy of the baryon with an er- 
ror of less than 1.5% even for full occupation (y — 1) and 
arbitrary 7. In the chiral limit, the error is below 0.3% 
for all values of v. Even the LO effective action yields 
results no worse than 10-20% for the binding energy. 

The most important result of this section is the ef- 
fective Lagrangian (f3"9"|) for the massive GN2 model. It 
contains the corresponding Lagrangian (|19p for the mass- 
less GN2 model as a special case (7 = 0). The applica- 
tions to homogeneous matter and baryons show that our 
general scheme is correct and underline the systematic 
character of the expansion. In the derivation, we only 
had to compute standard one-loop Feynman diagrams 
with increasing number of external lines. Nevertheless, 
we demonstrated that an excellent approximation to the 
baryon mass can be obtained even in regions of the pa- 
rameter space where we had no a-priori reason to assume 
that our truncation of the effective action was justified. 



III. NAMBU-JONA-LASINIO MODEL (NJL 2 ) 

The Lagrangian of the Gross-Neveu model with con- 
tinuous chiral symmetry (NJL2 model) is 

C = + y [ W) 2 + feVO 2 ] ■ (44) 

In this model, the Dirac sea has a more dramatic impact 
on the hadrons of the theory than in the case of the 



GN2 model. This is already clear from the existence of 
massless mesons and baryons, despite the fact that the 
elementary fermions acquire a dynamical mass [2j. As 
compared to the GN model, one has to expect additional 
complications when integrating out the Dirac sea. In the 
present section, we therefore restrict ourselves to the LO 
calculation of the effective action. It turns out that the 
LO calculation in the NJL2 model is about as complex 
as the NNLO calculation in the GN2 model described in 
Sect. II. 

The general procedure and the diagrams to be con- 
sidered stay the same as before, except that each vertex 
can now be either 1 (scalar) or i 7 5 (pseudoscalar). Some 
diagrams vanish due to parity selection rules. 

At one-loop order (tadpole), the scalar contribution 
yields an effective mass just like in the GN2 model. The 
pseudoscalar tadpole vanishes, so that the one-loop con- 
tribution to C e g is again a standard mass term, 



= -miptp. 



(45) 



At two-loop level, we have to consider that the two ver- 
tices in Fig. 5a can both be either scalar or pseudoscalar. 
The first case is identical to the GN2 model and, af- 
ter resumming the scalar bubbles, yields the first term 
in Eq. (|11[) (higher order terms are discarded since we 
now work to LO only) . The second case involves a pseu- 
doscalar vacuum polarization. For the sum of the " + " 
tadpole and this correction, we obtain to 0(fc 2 ) 



Ng 2 
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4m 2 / 



(46) 

This should be compared to the corresponding scalar re- 
sult, Eq. ([5]). Using the same arguments for resumming 
the inner "— " bubbles as in the scalar case, we get (de- 
noting the pseudoscalar coupling by G e ff) 



G 2 ff (fc) 
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Invoking the gap equation ([5]) now yields the result 



G 2 cS (k) 
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(47) 



(48) 



The pseudoscalar bubble sum has produced the pole of 
the massless "pion". Actually, the difference between 
(|48|) and the leading order result for the scalar coupling 



g 2 ff = n/N can easily be understood: As already noted 
by Gross and Neveu [l| , the coupling constants g^pp and 
g a pp are identical due to chiral symmetry. In Eq. (|12p . 
we have interpreted the scalar effective coupling constant 
in terms of a cr-exchange, replacing the propagator by a 
constant, namely its value at k 2 = 0. However, this can- 
not be done in the case of the 7r-exchange. The limit 
/c 2 — > is singular and we must keep the leading momen- 
tum dependence, 



1 



4m 2 7r 
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(49) 



thus reproducing Eq. (148)) . Clearly, the pole in the ef- 
fective pseudoscalar interaction has drastic consequences 
for the structure of the effective theory. Since the ex- 
pansion around k = is singular in the NJL2 model, 
we have to reconsider the ordering principle behind our 
approach. Keeping the momentum dependence from the 
pion pole in Eq. (f48|) . the total two- loop contribution to 
the effective Lagrangian assumes the non-local form 



2nm 2 1 

-^—ip^ip—ip^ip. (50) 



As before we shall treat 00 and □ as being of 0(e). 
We would like to test our effective action against the 
multi-fermion bound states of Shei, in analogy to our test 
against the DHN baryon in the GN2 case. From these 
known bound states we infer that 0i750 is of 0(e 3//2 ). 
Hence both terms in Eq. ([50]) are of the same 0(e 2 ), 
which defines our LO approximation. 

At three-loop level we must again consider the graph 
in Fig. 5b. If all three couplings are scalar, it is of 0(e 3 ) 
and therefore NLO. We neglect it in the present case. 
The only other non-vanishing contribution involves one 
scalar and two pseudoscalar couplings. We evaluate the 
corresponding Feynman graph and resum the scalar and 
pseudoscalar bubbles into effective couplings, generating 
two pion poles. The resulting effective Lagrangian reads 



C 



(3) 
off 



87r 2 m 3 - 11 
N2 ipi^ip—ipip—ipy^ip. (51) 



Notice that this is again of 0(e 2 ). Due to the inverse 
powers of momenta or □, unlike in the GN2 case, there 
is no suppression of the three-loop contribution as com- 
pared to the two-loop contribution by a power of e. Even 
worse, one can identify a whole class of higher terms with 
arbitrary many loops contributing to the same order of 
e. This calls once again for a resummation, namely [add 
up the 2nd term in Eq. (J5(| and Eq. (j5T|) ] 



27rm 2 7 If inm - 1 



0i7 5 



27rm 2 



-0i750 



1 



□ 



01750- 



(52) 



We have summed up a special class of higher order dia- 
grams, singled out by being of 0(e 2 ). They can be de- 
scribed as follows: Draw a string of pseudoscalar bubbles 
between two 0i750 (valence) condensates. Then attach 
at most one scalar tadpole to any of the pseudoscalar 
loops (on either side). All diagrams generated with this 
prescription are summed up. If there are n scalar con- 
densates, this must be matched by n+1 pion poles. It is 
easy to convince oneself that any other way of attaching 
tadpoles to the bubble graphs is punished by a higher 
order in e. As a result of this discussion, we replace the 
effective action ([50| induced by the two-loop graph with 
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FIG. 8: Dressed pion propagator (dashed line), connected to 
an arbitrary number of scalar densities 0^) via cr-exchanges 
(wiggly lines). See Eqs. ([52} and J53}. 



The massless pion propagator has been superseded by 
the propagator of the pion in a scalar background po- 
tential and sums up higher loop contributions of 0(e 2 ) 
in £ e ff as shown schematically in Fig. 8. Notice that the 
pion self-energy appearing in the denominator of Eq. (|53|) 
has a form reminiscent of the Gell-Mann, Oakes, Renner 
(GOR) relation [18| which relates the pion mass to the 
bare fcrmion mass and the chiral condensate, 



N 



<0V>. 



(54) 



However, in Eq. (|53p . m is the physical fcrmion mass, 
and ipip refers to the (x-dependent) valence contribution 
to the condensate only. Nevertheless, the formal similar- 
ity is striking and may point to a deeper physics reason 
behind our rather technical resummation. 

Due to the troublesome occurrence of inverse powers 
of e in the pseudoscalar effective coupling, we have to go 
to even higher order in the loop expansion in order to 
identify all LO terms. At four-loop order, we had to con- 
sider Fig. 5c for the GN2 model (where it was NNLO). If 
all four vertices are scalar (like in the GN2 model), it is 
indeed of 0(e ) and can be discarded for our present pur- 
pose. If two vertices are scalar and two are pseudoscalar, 
after resummation it will be of 0(e 3 ) and hence still neg- 
ligible. The only 0(e 2 ) contribution is the diagram in 
Fig. 5c with four i7s vertices. The four pion pole terms 
together with four pseudoscalar condensates conspire to 
give 0(e 2 ). The Feynman diagram calculation yields the 
effective action 
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According to our previous treatment of the two-loop di- 
agram, we should once again replace every denominator 
□ by □ — ^p0?/>, thereby summing up another infinite 
set of diagrams of 0(e 2 ) but containing multi-fermion in- 
teractions of arbitrary order, 



C 



(4) 



327T 3 TO e 
TV 3 



1 



□ 



01750 



(56) 



Finally we turn to the five-loop graph in Fig. 5d, which 
contributed in NNLO in the GN2 case. The interesting 
case is the one where there are four pseudoscalar ver- 
tices and one scalar vertex, the latter connecting the two 
" — " bubbles. The calculation including all necessary 



b) 



b) 




FIG. 9: a) Effective direct 4n interaction corresponding to 
Eq. (|56[) . b) Effective 4ir interaction via er-exchange leading 
to Eq. I|57p. The two diagrams cancel exactly due to the van- 
ishing of the irn interaction at zero momentum, in accordance 
with low energy theorems. 



FIG. 10: a) Effective fermion-fermion interaction through 
dressed 7r-exchange, giving rise to the non-local term in C e s, 
Eq. (|58p. b) Effective fermion-fermion interaction through 
(T-exchange, yielding the local {Tpip) 2 -terra in C c g. 



resummations yields 
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It yields exactly the same magnitude as (|56|) but the op- 
posite sign, so that the four- and five-loop terms cancel 
in the NJL2 model. In our opinion this cancellation is an 
expression of the low energy theorem according to which 
the 7T7T interaction at zero momentum should vanish . 
Here this comes about as a result of a cancellation be- 
tween the 47T vertex and the process 7T7t — * a — > tttt in- 
volving the 7T7TCT vertex, see Fig. 9. It is interesting that 
the NJL2 model respects the low energy theorems in the 
large N limit, in spite of the usual strong reservations 
against Goldstone bosons in two dimensions [20| . 

We have now identified all terms of 0(e 2 ) in the effec- 
tive Lagrangian. As announced, the calculation was more 
involved than in the GN2 model. Whereas in the GN2 
case we had to deal with four-, six- and eight-fermion in- 
teractions in LO, NLO and NNLO, here we had to sum 
up terms corresponding to 2n-fermion interactions with 
arbitrary n already in LO. The reason is the massless pion 
pole which induces inverse powers of e. It is responsible 
for the non-locality of the LO effective interaction for the 
NJL2 model which nevertheless has a rather simple final 
form, 
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The pseudoscalar term is associated with the fermion- 
fermion interaction through 7r-exchange, see Fig. 10a. 
The pion propagator has a long range and is dressed as in 
Fig. 8. The scalar term arises from cr-exchange, Fig. 10b, 
and is of zero range in LO. It is common to the GN2 and 
NJL2 models. 

At first glance, the non-locality of £ e ff is irritating, 
since it is expected to make the solution of the Eulcr- 
Lagrange equation much harder. However in the present 
case, the particular structure of the non-locality allows us 
to trade the non-local effective Lagrangian for a local one, 
at the cost of introducing an elementary pseudoscalar 
field U(x). The local effective Lagrangian equivalent to 



58]) is given by 



eff 
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Since we are in the large N limit, n can be treated as 
a classical field. The equation of motion for n following 
from ASS]) is 
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(60) 



Upon solving this inhomogeneous equation formally and 
plugging the result into Eq. (|59p we indeed recover the 
non-local effective Lagrangian (|58[) . So we wind up with 
a local field theory containing both fermions and an ele- 
mentary 7r-meson, the latter described by the field H 

Since the derivation was quite complicated (in particu- 
lar identifying all LO terms), it is again crucial to test C c s 
against exact results for the NJL2 model. For this pur- 
pose, we choose the multi-fermion bound states of Shei 
[6j. They are analogous to the DHN baryon of the GN2 
model in the way they were derived (inverse scattering 
theory) , but in our opinion carry vanishing baryon num- 
ber and should be regarded as "baryonium" states (Rl2l|. 

Our strategy is as follows. Take the valence spinor 
wave function of Shei in the notation of Ref . 01 , 



m[cos_g| 1 / cos 6>/2 \ 
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with 
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and energy Eq = —msm8. This yields the (valence) 
condensates 
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Eq. ([rjO]) for n becomes 
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(64) 
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and, to leading order in v, is solved by 



n = - 



1 + e 2 « ' 



(65) 



The Euler-Lagrange equation for tpo is equivalent to the 
"no-sea" Dirac-HF equation, 



(66) 



with scalar and pseudoscalar potentials 



S(x) 



-IT = mil 



2n 2 v 2 
1 + e 2 « 



4ir 2wv 

Pix) = —\ — mil = m TTT- 

V ' V N 1 + c 2 « 



(67) 



We have kept the LO terms only. We find that ipo sat- 
isfies the HF equation up to corrections of 0(^ 3 ). The 
Hamiltonian density is given by 
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Taking its expectation value in the valence state, we in- 
sert the above results for ipo and II and integrate over x. 
The derivative term gives zero, the mass term has to be 
treated in NLO in v and yields Nmv — Nmir 2 is 3 /2, the 
four remaining terms give Nmir 2 v 3 /3 to LO. The total 
result for the mass adds up to 
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dxH = Nmv 1 



2 2 
7T V 



(69) 



in agreement with the first two terms of the Taylor ex- 
pansion in v of the exact result, 



Nm 

M = sin irv. 



(70) 



Since we cannot expect more from a LO calculation, the 
test was successful. 

Now consider the question of baryon number. One 
interesting aspect of the NJL2 model is the fact that a 
topologically non-trivial mean field can induce fermion 
number in the Dirac sea [HI, l23l | . Since this is a physi- 
cal effect but we have eliminated the Dirac sea, it must 
show up somewhere else in the effective theory. To un- 
derstand what is happening, we treat the fermion density 
tp'tp — r ip^ ip as an effective operator, in analogy with 
the scalar or pseudoscalar densities above. To this end 
we write down the dressed tadpole with a 7 vertex and 
again resolve it into positive and negative energy contri- 
butions. The diagrams involving only " — " terms add 
up to the divergent density of the (vacuum) Dirac sea 
and can be subtracted. The " + " tadpole represents the 
density contribution of the valence level, p va \ = nipQipo, 



FIG. 11: Illustration of the effective fermion density, Eq. (|71|) . 
a) Valence contribution, b) Vanishing sea contribution due to 
the scalar interaction, c) Induced fermion density due to the 
pseudoscalar interaction. 



see Fig. 11a. Its first correction is driven by the two- 
loop graphs in Fig. 11. If the upper vertex is 1 (scalar 
" + " loop, Fig. lib), the vacuum polarization bubble is 
the time component of a four-vector and therefore pro- 
portional to k° (where k is the external momentum). 
Since the self-energy insertion is static, this term van- 
ishes. This is the reason why we did not have to worry 
about induced fermion number in the GN2 model. If the 
upper vertex is 175 (pseudoscalar " + " loop, Fig. 11c), 
the " — " bubble behaves like the space component of 
a four- vector and is proportional to k . We can replace 
the pseudoscalar " + " bubble by the full pseudoscalar 
self-energy P, thereby summing up all the LO diagrams 
as explained above. The result for the induced fermion 
density in coordinate space then assumes the simple form 
Pind = \J N /ird x Tl, where the derivative reflects the above 
mentioned /^-dependence. The total fermion density in 
the effective theory consists of the valence fermion den- 
sity and the induced one, 



p = Pvai + Pind = n%ipo + 



a*n. 



(71) 



This allows us to correctly compute the fermion density 
without Dirac sea. If we insert the expressions for the 
Shei bound state, we get an exact cancellation between 
the two contributions, in agreement with two indepen- 
dent recent results [1, [H|. 

Let us now switch on the bare fermion mass in the 
NJL2 model. As discussed above, this modifies the vac- 
uum gap equation, see Eq. (pJB]) , and thereby the effective 
couplings. For the scalar coupling the corresponding LO 
result can be read off from the (4>ip) 2 term in Eq. 
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An analogous calculation for the pseudoscalar effective 
coupling yields 



„ = 4m 2 7T 
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The massless pion pole gets replaced by a massive one 
with the pion mass 



m 2 = 4m 2 7, 



(74) 



in agreement with the leading order pion mass prediction 
near the chiral limit (GOR relation 0, [l8j]). The next 
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steps depend on the regime one is interested in. Two 
cases are particularly simple and will be considered here: 
i) m 2 and k 2 are both of 0(e), and the scalar and pseu- 
doscalar condensates are 0(e) and 0(e 3 / 2 ) respectively as 
for the Shei bound state, ii) 7 » 1, i.e., ir and a masses 
are comparable and of 0(m), cf. the exact relationship 
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4m 2 



(75) 



In case i) we can neglect 7 in the scalar effective cou- 
pling multiplying (V>?A) 2 since it is of higher order than 
0(e 2 ). The same diagrams contribute to the pseudoscalar 
terms as before, except that the pion proapagator be- 
comes massive everywhere, □ — > □ + m 2 . The cancella- 
tion in the irir scattering amplitude illustrated in Fig. 9 
is upset by the bare mass term, as expected from chiral 
perturbation theory and low energy theorems. However, 
since the correction gets an additional factor of 7, these 
terms can again be ignored to LO. As a result, the non- 
local effective action (|55|) of the massless NJL2 model 
now gets replaced by 
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with m 2 defined in Eq. (I74| . This can again be converted 
into a local effective action by introducing an elementary 
pseudoscalar field II. The only modification as compared 
to Eq. ([59)) is the appearance of a mass term for II, 

4° ff c = ^- m )^ + ^(^) 2 + i^n^n (77) 

27TTO - , „n 1 o „n 

ibibtt 2 - -mitt 2 . 

N 2 ff 

This opens up the possibility to extend Shei's bound state 
to finite bare fermion masses. Of course one would then 
have to verify that our assumptions concerning the order 
of e of the valence condensates are fulfilled. 

In case ii) we go into a regime where both tt and a 
propagators may be treated as point-like. The scalar and 
pseudoscalar effective couplings become equal (to LO), 



.9e ff = G 2 ff = — , 7 »1. 



(78) 



All higher order corrections are suppressed by inverse 
powers of 7, and the LO effective Lagrangian assumes 
the same structure as the original Lagrangian (apart from 
the mass term and the replacement of g 2 by the effective 
coupling) , 

£ off = i> (i? - m) ij> + ^- [W) 2 + i^l^f] • (79) 

In this regime, the LO effective Lagrangian of the GN2 
model reduces to 



ip (i^ — m) ip ■ 
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(80) 



These last two equations may be considered as the ex- 
treme heavy-fermion (or non-relativistic) limit of the 
original field theories. 

Finally, we comment on some early work which is rel- 
evant in this context. In the 70 's, there was some debate 
about the use of "classical calculations" for fermions in 
field theory. These were mean field calculations without 
Dirac sea, solving non-linear Dirac-HF equations [[J, [l(| 
exactly as in our "no-sea" effective theory. What was 
done at that time simply amounted to use the original 
Lagrangian with a mass term added and replacing the 
bare coupling constant by an effective one. Shei [f| noted 
that this type of calculation gave reasonable results for 
the GN2 model, but was completely off in the case of 
the NJL2 model. We can now understand the reasons 
behind this observation, using the framework of effec- 
tive field theory. In the GN2 model, the lowest order 
effective Lagrangian happens to agree with the ad-hoc 
prescription used in Refs. [^, [h|. This explains why 
it gave useful results, at least for small filling fraction. 
Incidentally, another example of a "no-sea" calculation 
for this model can be found in Ref. Q where the effec- 
tive coupling constant g 2 s = ir/N was determined "phe- 
nomenologically" , see Eq. (4.17) of that paper. We have 
now derived this number microscopically and are able to 
improve the calculation in a systematic fashion. In the 
NJL2 model, the naive recipe for writing down the ef- 
fective action breaks down due to the pion pole. One 
cannot ignore the singular momentum dependence of the 
effective pseudoscalar coupling (~ 1/fc 2 ) which gives rise 
to long range non-localities. This is the reason for the 
observed discrepancy between short range "no-sea" and 
long range semi-classical mean fields. It is amusing that 
in the heavy-fermion limit, we do recover an effective ac- 
tion which agrees with the simple recipe of adding a mass 
term and introducing an effective coupling, see Eq. (|T9|) . 
Hence the technically rather nice classical calculations 
done in the 70's for the chiral limit of the NJL2 model 
(where they had to fail, as we now understand) can be 
"recycled" to solve the NJL2 model in the heavy-fermion 
limit 7 > 1. In this regime, the approximations can be 
justified by effective field theory. In the representation 
of the 7- matrices (|33)l . the bound state found by Lee et 
al. fiol ] has the spinor wave function (filling fraction v 1 
energy E ) 



ipo 



(m - E ) 



2G 



cosh 2 £ + a 2 sinh 2 £ 



a sinh £ + cosh £ 
a sinh £ — cosh £ 



G = — , a = tan—, 
Eq = mcosG, £ = mxsinG 



(81) 



We have checked that it satisfies the (classical) non-linear 
Dirac equation obtained from the effective Lagrangian 
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(|79[) exactly. The baryon mass is given by 

sinG / -k 2 v 2 \ 

M B =nm-^-=nmll-^ T + ...J (82) 

Due to the truncation of the effective action, only the two 
listed terms of the series expansion can be trusted. At the 
same time the induced fermion density gets suppressed 
by a factor of 1/7 as compared to the valence density. 
The corresponding result for the GN2 model with effec- 
tive action ([80]) can be read off Eq. (J43J) (for 7 > 1) and 
agrees with (|82|) to this order in v. As a matter of fact, 
the spinors V'o also become equal in the limit 7 — > 00, 
as can be seen by expanding in powers of v. The reason 
can be traced back to the fact that the ratio of pseu- 
doscalar to scalar condensates decreases like I/7. Thus 
we learn that in the heavy- fermion limit, the difference 
between discrete and continuous chiral symmetry of the 
four- fermion interaction becomes less and less important, 
and baryons of the GN2 and NJL2 model approach each 
other. As far as the NJL2 model is concerned, there is 
so far no independent information available about this 
regime these findings could be compared to. 

IV. SUMMARY AND CONCLUSIONS 

Using exactly solvable four-fermion models as theoret- 
ical laboratory, we have demonstrated that "integrating 
out the Dirac sea" is indeed a viable concept. Tractable 
effective Lagrangians could be derived, valid for multi- 
fermion bound states with weakly occupied valence level 
or large bare fermion masses. Whereas the original 
semi-classical calculations of bound states in GN2 and 
NJL2 models required sophisticated and highly special- 
ized techniques to deal with the Dirac sea, the effective 
theories amount to solving the classical Euler-Lagrange 
equation, here a non-linear Dirac equation. One might 
suspect that the difficulties of the full non-perturbative 
calculations come back in the derivation of the effec- 
tive Lagrangian, but this is not quite so. Although 
this derivation has non-perturbative features in the form 
of unavoidable resummations, the actual calculation is 
based on standard one-loop Feynman diagrams and thus 



perturbative. This suggests that similar techniques can 
also be applied to more realistic, higher dimensional the- 
ories. 

In the GN2 case, we were able to derive the NNLO ap- 
proximation in some small expansion parameter. Since 
the exact DHN baryons and their generalization to fi- 
nite bare fermion mass are known, we could show that 
the resulting effective theory works quantitatively. Thus 
baryon binding energies can be predicted "classically" in 
the full range of filling fraction and bare mass at the 
1% level, quite unexpectedly for us. As far as the for- 
malism is concerned, the main lesson we learned is that 
the bare parameters of the original theory (mg^g 2 ) dis- 
appear owing to resummations. It would not make sense 
to truncate the procedure at any fixed number of loops. 
All questions of renormalization and UV divergences are 
then properly dealt with and do not show up any more 
in the effective theory, which is of course finite. 

In the NJL2 case, when applying the same strategy 
we ran into a new difficulty. The pseudoscalar effec- 
tive coupling develops a singularity at k 2 = due to 
the massless pion pole. Although one can again identify 
the LO of a systematic expansion in a small parameter, 
this necessarily causes a long range non-locality of the 
effective Lagrangian. The purely fermionic, non-local ef- 
fective theory can be cast into a more convenient "almost 
local" form by introducing an additional elementary pion 
field. The resulting theory contains elementary fermions 
and pions but is fully consistent, judging from the deriva- 
tion and the comparison with Shei's bound states. The 
role of the massless pion explains the conspicuous differ- 
ence between short range mean fields in the GN2 model 
and long range mean fields in the NJL2 model, which 
had caused headaches in early, more naive attempts to 
treat fermions classically. Moreover, we could show how 
the important phenomenon of "induced fermion number" 
can be recovered in a "no-sea" effective theory. 

Finally, notice that no attempt was made to reduce the 
Dirac equation to a non-relativistic Schrddinger equation. 
Although this would be possible, we feel that it would 
only make our formulae more messy. Our main focus here 
was not on relativistic kinematics, but on the dynamical 
effects of the Dirac sea. 
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